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BOUNDARY VALUE PROBLEM FOR TWO-LAYER 
HALF-PLANE 


O. YAREMKO 


ABSTRACT. In the article the formulas for the modeling of conservative fields in 
piecewise infinite plate with a thin inclusion found. The accuracy of the found 
formulas is of order equal to the thickness of the outer layer. The problem for 
higher accuracy solved by means of asymptotic formulas. 


1. INTRODUCTION 


Problems about structure of the conservative field of a two-layer flat plate leads to 
a separate system of the Laplace’s equation with a boundary condition to Dirichlet 
and internal coupling conditions [1], [2]. In the work the method of transformation 
operators is applied. Yaremko O. E. developed this method [1], [2], [3], [4], [5]. As 
a result the conservative field for the two-layer piecewise plate can be interpreted as 
a deformation of the conservative field of a homogeneous piecewise plate. Thus the 
deforming transformation operator is written out. The deforming transformation 
operator is recursive, and convenient for practical realization on the computer. 
In this case, use the analytic representation of a field by means of the deforming 
transformation operator is difficult. 

We consider the case of small thickness | cross-border layer of the plate. Let’s 
consider that the physical properties of the layers very different, that is the ratio 
of the thermal conductivity of the layers close to zero. 

The method of the deforming transformation operator gives algorithm inconve- 
nient for realization. Formulas which have no these defect are found based on the 
Euler—Maclaurin formula. 

Numerical methods of research lead to ill-conditioned [8] systems of the linear 
algebraic equations. 


2. THE PLANE CASE 


2.1. The half-plane with internal coupling condition. Modelling of potential 
fields for a semi-limited infinite plate with thin inclusion leads to the Dirichlet prob- 
lem for two-layer half-plane with coupling conditions on the direct. Let’s consider 
a boundary problem about the solving of separate system of the Laplace’s equation 
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uz 3u 
(2) a2 ggr + gyr = 9: Eyle =o <y <o 
by the boundary condition for direct x = l : 


(3) ui(z,y) = û(z, y), £ = 0 

by the internal coupling conditions [1], [2], [3], [4], [5] for direct x = l : 
ð ð 

(4) u(x, y) = uaz, y), k= (z, y) = S (2,4), k > 0, 
Ox Ox 


Here function & = t(x,y)- is harmonic in the right half-plane, H = {(z,y) : 0 < 
x,y E€ R} and continuous in {(x,y):0 < x,y € R} and the conditions are met 


û = 0(y z? + y?), 


û(0, y) 
d 
f Tae yY <0 





—oCo 
The problem (1)-(4) is ill-posed [8] in the case of a thin inclusion, that is, small 
value 1 , and in the case of a strong heterogeneity, that is, large differences in the 
parameters a,,a2. The solution of this problem by known methods is accompa- 
nied by considerable computing difficulties. New analytical methods to solve this 
problem it is necessary to develop. 


2.2. Thin strip. Function @ = éi(x,y) is a harmonic function in the right half- 
plane H = {(x,y) : 0 < x,y E€ R} continuous in {(2,y) : 0 < x,y € R} and the 
conditions are met 





Let’s consider the Dirichlet problem for Laplace’s equation in the strip Hı = {(2, y) : 
0<a<lyeER} 


8u 8u 
©) Oe t OP 
with boundary conditions 
(6) u(0, y) = a0, y), u(l, y) = 0. 
The problem (5)-(6) is ill-posed [8] at small value 1. 


= 0, (x, y);0 < xz < l,—00 < y < œ 


3. THE AXIAL CASE 


3.1. The circle with the internal coupling conditions. Let the function û = 
a(x, y)- is a harmonic function in the unit circle, B = {(x,y) : £? +y? < 1} and 
continuous in its closure B. The boundary problem about the solving a separate 
system of the Laplace’s equation consider in the circle B 


Au = 0, (x,y) € Kr, 


(7) Aue = 0, (x,y) € Br 
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Here Kp - annulus with the smaller radius R and larger radius 1, Br - the circle 
of radius R. Let the boundary condition on the circle satisfied S: 
(8) u(x, y) = û(x, y), (a, y) € S 
internal coupling condition on a circle of radius R: 


u(x, y) = u2(£,y), 


(9) kLoui(z, y) = Lou2(x, y), k> 0, (x,y) = SR, 
o o 
Lo = £— —. 
0 = Tag TA D 


Let’s notice that the problem (7)-(9) is ill-posed [8] at values R ~ 1. 


3.2. Thin annulus. Let the function û = û(x,y) - is a harmonic function in the 
unit circle, B = {(x,y) : 2? +y? < 1} and continuous in its closure B. In annulus 
Kpr- with the smaller radius R and larger radius 1, consider the Dirichlet problem 
about the solution of the Laplace’s equation 


(10) Au; = 0, (x,y) € Kr, 
by the boundary conditions on the circumferences S, SR: 
(11) u = û, (x,y) E€ S; u = 0, (x,y) € Sp. 


The problem (10) is ill-posed [8] by the values R = 1. 


4. TRANSFORMATION OPERATORS 


In the works of the author [1], [2], [3], [4] and [5] the concept of operator trans- 
formation J, J: û — u, allowing a known solution ù of the model problem to solve 
u any of the four above-listed problems determined. 

The half-plane with the internal coupling conditions 

Transformation operator J has the form 


J: û > u,u = x(Hı)uı + x(H2)u2, 








where: 
Leg 1—-k 
je; ——$———— 7 f ) 7 j l 
u(y) la) (ate +2449) - Rae- eai), O<e<t, 
Die Leo Pag 
= 7 j l 
(12) u(x, y) a) «(Ste )+1+ 2y), <Q, 
Ay a2 
k= —— 
2 a4’ 


x— is the characteristic function of the set [6]. 
Thin strip 
Transformation operator J : ti —> u, can be determined by the formula 
(13) u(a,y) = >> (ala + 2lj, y) — û(2l — x + 2lj,y)); O<a<l. 
j=0 





Circle with the internal coupling conditions 
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Transformation operator is expressed by the formula: 


J : û —> u,u = x(Kr)uı + x(Brlu2, 








where: 
2 /1-k\f ; ~ l-k, (RIĦ R2+? 

14 = a 2lj 2lj a 
(14) u(z,y) 2 (=) (aur yR") TERE = y ; 

Ik EAT LEEN sa mee 

û J 2j 
waleo) = FFT! Aa a(R, yR”), 

x — is the characteristic function 7 the set [6]. 


Thin annulus 
Transformation operator J : û > u, can be determined by the formula 


Co 


(15) u(a,y) = >> (« (2R, yR) —u = ZN, 


j=0 Y 





Formulas (12) and (14) are useful only if the values k ~ 1. In the case of the two- 
layer body, which components are sharply different properties, that is at values k 
close to zero, and for large values k, transformation operators cannot be used, in 
view of the slow convergence of the series at I” xl. 

Formulas(13) and (15) cannot be used in view of the slow convergence of the 
series at 1 ~ 1 and R ~ 1, respectively. 

The purpose of our research is to consists in creation of a design of transformation 
operators comfortable in these cases. 


5. AUXILIARY RESULT 


Lemma 1. If the function y = f(x) is defined on the line segment [0,1] and has 
limited variation [10] on the line segment [0,1], then for each value R the following 
assessment for the difference of integral and its integral sum is carried 


1 


< 1 
[to-ny mL R| < m VC) 
J : 


Proof. Let’s spread integral into the sum (n +1) of the summand: 
1 n-i 1 
fa flejde = X f £t f(R” x)dr + fx! f(R”x)dxrVariable substitution in 
0 j=0 R2 0 

the last integral results in equality: 


R?” 


/ a f(a m priya + f tod 


I= OR2 0 


; 1 
Considering that |f (Rx) — f(R%)| < V(f), We obtain the required assessment: 
0 
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) dz 1 
<È f pria- re 7 


j= Op2 


<V(f)In 


1 
V(f) — variation of function [10] f on the line segment [0,1]. The flat analog of 
0 


the Lemma is well known [8]. 

Lemma 2. If the function y = f(x) is defined on the ray [0;co) and has limited 
variation [10] on the ray [0; co); then for each value 1 the following assessment for 
the difference of integral and its integral sum is carried 


[far a] < UF) 
0 j=0 


The Bernoulli numbers [11] are defined using the generating function 


The Euler—Maclaurin formula applied to the function f(2lx) on the ray [0;co) results 
in equality 





oo Pa œ 2k-1 
Dren = g f toat LO — SPB peng, 
0 





j=0 k=1 

3 17 1) & (2)*- a 
Leys Af paar E EGG) 
j=0 In z 3 k=1 


Corollary 1. When k € (0,1) completed asymptotically 


°° CENI 
Shia) f(a + 2j) = 
j=0 





1 
2l 


i 





fe pee + eae-+ £2) 5 O Bu paa pia) 
0 k=1 


When k € (1,co)completed asymptotically 








Sf LRN why, = o 2 
(Fz) tora at (g?* — 1) EFF 4 F(a), 
j=0 k=1 
The operator Ləp has the form 
d 
Lop, = 2h + —. 
2} t 


Corollary 2. When k € (0,1) completed asymptotically 


oO a j i 
(a) for 
0 


j= 





1 1 ‘eae 


= = perce: a H >D R? hl Lahey), 


~ in & 
R 0 








When k € (1,00) completed asymptotically 


© = j ; r (In 4,)2*-1B, r ai 
S (TA) tern AD y Te ok i e). 


j=0 





The operator Ləp has the form 


d 
Lon = 2h + —. 
2} Ta 


6. MAIN RESULT 


6.1. The plane case. Study the case of a thin shell, that is the case of a small 
thickness of the final layer 1. Physically this means that the coefficients of heat 
capacity layers are very different. 
Theorem 1. Let the condition is fulfilled 

l-k 2hl 

(ee e? 
and the value k is small. For component u; and uz of the problem solution (1)-(4) 
approximate formulas are valid 


1 — e2hl 
Uz X ay ts» 
1 e2h 
u, ~ Sts (wy) — zr ĉel- x,y), (2, y) € H2, 


Where ûz - solution of the problem of Robin [12] with the boundary condition 


ð 
2hi3 + —tiz = ù, (x,y), £ = 0. 





On 
The assessment is fair 
1 —e2nl F zs 
A — | alet 2,y)de— us(e.y)] < (0-2) Peale + ey) 
0 


0 


Proof. Let’s apply the lemma 2. As a result we obtain the formula: 


2k 1 ‘ner 2k cha 

eae a = Zo £ : 

aan le (a + €,y)de — u2(x,y)| < Rei ¥’ (a +e,y)) 
0 


We will express k through h. Get the assessment, which proved for u2. The formula 
links solutions the Dirichlet problem and the Robin problem from [12] 


Co 


uz = — / ealz +e, y)de. 
0 
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This theorem allows the physical interpretation: component of the solution uz with 
the accuracy to a numerical multiplier it is approximately equal to the solution of 


the third homogeneous boundary problem with boundary condition 


ð 
2hûz + —tiz = ti, (x,y), £ = 0. 


On 
Similarly 
Theorem 2. Let the condition is fulfilled 
l-k _ aahh 
p 


and the value k>>1. For the problem solution (1)-(4) approximate formulas are 
valid 
X {= e2hl 1 7 ao J 
ui X T+ eh ((ai3(2, y) — "ig (a2 + ,Y)) + tis) , 


(z + 7"! (ûa (2l — z, y) — e” ûs (4l — z, y))) 7 


S j= e2P! 1 7 Pic j 
Ug ay (ua(eu) —€ ti3(a + ,y))- 
Theorem 2 is proved on the basis of the Lemma 2. 

Thin strip 

We apply Lemma 2. For the solution of the Dirichlet problem in the strip (5)-(6) 
assessment is fair 


e bevy) -êl -r +e) 
0 











e—ulz,y)| <V, 
0 


2l 
here 


- variation of function [10] taken from a variable € on the interval [0;co). Thus, the 
following theorem is obtained 

Theorem 3. The solution of the Dirichlet problem in the strip (5)-(6) in the 
strip can be found on the approximate formula 


û2(x, y) = ta (21 =T, y) 
2l f 
Where ûz - the solution of the Neumann problem for the equation Laplace’s in a 
circle with the boundary condition 
ûz 


Ox 





ula, y) ~ 


=t,x =0. 


6.2. The axial case. The circle with with the conditions internal coupling. Study 
the case of a thin shell, that is the case of a small thickness of the outer layer 1-R. 
Let’s consider also magnitude k of the small, that is consider the case when the 
coefficients of heat capacity layers differ greatly. Assuming 


1-k 
— = R™. 
1+k 


Formula follow from Lemma 1 
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2k 1 2k 1 
G z J ales, ye)de — u2(x,y)| < ERN y (e"û(xe,ye)), 
0 
1 
1— R?” h Lah h 
i fe ~"&(ae, ye)de — uz(x,y)| < V (e"û(ze, ye)) (1 — R”). 
In RZ n 0 


Therefore, we come to the next result. 
Theorem 4. For the solution of the problem (7)-(9) u1, u2 approximate formulas 
are valid 











ui X ’ 


zT y 
This theorem allows the physical interpretation: component of the solution uz with 
the accuracy to a numerical multiplier it is approximately equal to the solution of 
the third homogeneous boundary problem with boundary condition 


1. R? — /R R 1- R?” 
In I ûz(x, y) In 3 U2 X u3. 


1 1 
RZ R In Re 


2hiiz + 2 i = ù, (x,y) E€ S. 
It is necessary to consider the formula from [3], The formula links solutions of the 
first and third homogeneous boundary problems: 
1 
a3 = [et tiles, ye)ae. 
0 


Study the case of small values of thickness of the external layer 1 — R ~ 0, thus we 
consider coefficient k>>1. If the expression component us rewritten in the form: 








Qo a fk-1\" as as 
wl) = E> (FS) i(cR* ,yR*)— 
j=0 
2k pol /k-1\”. ; 
Saag rer (7) i(aR? RY ,yR? RY”), 
j=0 
define the number h using the formula 
k-1_ 2h 
1+k 


and lemma 1. We get the approximate formulas for component of the problem 
solution (7)-(9): 


1-R* 1 
1+ R? 21n gs 


2 p2 4 p4 
Het a o o 
T y T y 





((a3(x, y) — R? ûs (R?a, R?y))) + 


Uy ~ 





1- R” Qhe (p2, p2 
US ani (ûs(x, y) — R” ûs (Rx, R?y)). 
n Rz 


Thin ring 
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Let - the function & = û(x,y) - is harmonic in the unit circle, B = {(x,y) : 
x? +y? < 1} continuous in its closure B. In the ring Kp- with the smaller radius 
R and larger radius 1, consider the Dirichlet problem (10)-(11). 

We apply Lemma 1, we obtain the following assessment of solutions of the Dirich- 
let problem in the annulus (10)-(11) 





here 





1 1 . . R? R? 
V = V | û(xe,ye)— û E, —E 
0 0 “uy 


variation of function [10] taken from a variable ” on the line segment [0; 1]. Thus, 
the following theorem is valid. 

Theorem 5. For the solutions of the Dirichlet problem in the annulus (10)-(11) 
approximate formula is valid 





ula, y) = 


Where tz - the solution of the Neumann problem for the equation Laplace’s in a 
circle with the boundary condition 
Ou 


an = ti, (x,y) € S. 


7. CONCLUSION 


From the results given in article follows, that the accuracy of the found formula 
is order equal to the thickness of the outer layer. Therefore, obtaining the formulas 
of higher accuracy class represents an important interest. This problem is solved by 
using asymptotic formulas presented in corollary 1, 2. The problem of distribution 
of the results on the multilayered plates and also on the boundary problems with 
the more general boundary conditions is actual. 
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